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,i— > Abstract 

5/3 The paper studies the interaction of a longitudinal wave with trans- 



verse waves in general isotropic and unconstrained hyperelastic ma- 
terials, including the possibility of dissipation. The dissipative term 
'j^ chosen is similar to the classical stress tensor describing a Stokesian 

'^ fluid and is commonly used in nonlinear acoustics. The aim of this 

5 research is to derive the corresponding general equations of motion, 

Q valid for any possible form of the strain energy function and to in- 

vestigate the possibility of obtaining some general and exact solutions 
CN to these equations by reducing them to a set of ordinary differential 

^ equations. Then the reductions can lead to some exact closed-form 

\^ solutions for special classes of materials (here the examples of the 

CN Hadamard, Blatz-Ko, and power-law strain energy densities are con- 

^ sidered, as well as fourth-order elasticity). The solutions derived are 

^T in a time/space separable form and may be interpreted as general- 

~^ ized oscillatory shearing motions and generalized sinusoidal standing 

T-H waves. By means of standard methods of dynamical systems theory, 

^ some peculiar properties of waves propagating in compressible materi- 

als are uncovered, such as for example, the emergence of destabilizing 
effects. These latter features exist for highly nonlinear strain energy 
^ functions such as the relatively simple power-law strain energy, but 

they cannot exist in the framework of fourth-order elasticity. 



1 Introduction 

Nonlinear elastic wave propagation is a subject of considerable interest for 
many natural and industrial applications such as seismology [S3], soft tissue 
acoustics [Sn|, or the dynamics of elastomers ^Tj. A detailed study of the 
theoretical issues associated with wave phenomena forms the basis of our 
understanding of important non-destructive and non-invasive techniques of 
investigation such as for instance, the technique of transient elastography for 
the analysis of soft solids [20] . The mathematics of wave phenomena is still 
an active subject of research where many outstanding problems are waiting 
for a definitive systematic treatment. Many papers are devoted to the study 
of nonlinear elastic wave propagation and a recent summary on the status 
of contemporary research on the subject can be found in the authoritative 
review by Norris [U], while an account of the mathematics of hyperbolic 
conservation laws can be found in a book by Dafermos [23] • 

It is important to note that in most studies of nonlinear acoustics, the 
investigation is usually restricted to weakly nonlinear waves. Indeed, for 
materials such as metals, rocks, or ceramics, the ratio of the dynamic dis- 
placement to the wavelength is a small parameter and so, the theories of 
third-order elasticity or of fourth-order elasticity seem to be sufficient to ac- 
count for the nonlinear effects observed in experimental tests. The situation 
is however completely different when we consider elastomeric materials such 
as those used for vibration isolators or automobile tires, and when we con- 
sider biological materials such as arterial walls or glands, under physiological 
or pathological conditions. The deformation rates occurring in these mate- 
rials are so large that their behavior will differ significantly not only from 
the behavior predicted by the linear theory of elasticity, but also from the 
behavior predicted by the weakly nonlinear theories. It should therefore be 
important to be able to derive and to investigate the general equations gov- 
erning the propagation of finite amplitude longitudinal and transverse waves 
in the context of full nonlinear elasticity, which recovers the weakly nonlin- 
ear case as a special case. Moreover, recent researches in the constitutive 
behavior of rubber-like material indicate that the stiffening effect, a peculiar 
but real phenomenon occurring at very large strains, cannot be described 
accurately when the full polynomial strain energy function is approximated; 
also, high order polynomial theories may give rise to numerical difficulties 
in the fitting of the material parameters with the experimental data and to 
artificial instability phenomena (see Pucci and Saccomandi [5Tj and Ogden 



et al. 

These considerations have led to the present paper which is devoted to 
the study of the interaction of a longitudinal wave with transverse waves 
in general isotropic and unconstrained hyperelastic materials and this also 
when dissipation is taken into account. The aim of our research is first, to 
derive the corresponding general equations of motion, valid for any possible 
form of the strain energy function; second, to investigate the possibility of 
obtaining some general and exact solutions to these equations by reducing 
the specific field equations to a set of ordinary differential equations; and 
third, to provide, using these reductions, some exact closed-form solutions 
for special classes of materials. 

To the best of our knowledge, previous studies taking systematically into 
account the full nonlinear equations of motion have generally been restricted 
to the formal theory of singular surfaces and of acceleration waves. These 
studies stem from the fundamental researches initiated by Hadamard [30] in 
1903 and were advanced mainly by Ericksen, by Thomas, and by Truesdell in 
the 1950s and 1960s (see for example Truesdell [58] and the review by Chen 
[22j|). They show that the conjunction of the nonlinearity in the material 
response and of the hyperbolic nature of the governing equations (in the 
purely elastic case) leads inevitably to shock formation. The results obtained 
are general and exact but they are formal results: in effect, no solution is 
found explicitly but rather, conditions are given that have to be fulfilled by 
solutions, if they exist. 

In contrast, the present article focuses on the explicit determination of 
solutions. Our investigation is directly related to the celebrated finite am- 
plitude elastic motions discovered by Carroll [HI [IHl [IB [121 II31 [II [IS] ■ His 
exact solutions are versatile in their fields of application because they are 
valid not only for nonlinearly elastic solids, but also for general viscoelastic 
solids y^, Reiner- Rivlin fluids [Ill[l5], Stokesian fluids tl4j, Rivlin- Ericksen 
fluids [I5], liquid crystals [I6], dielectrics [17], magnetic materials [H], etc. 
They also come in a great variety of forms, as circularly-polarized harmonic 
progressive waves, as motions with sinusoidal time dependence, as motions 
with sinusoidal space dependence, etc. Recently the present authors [25] ex- 
tended Carroll's solutions to the case of an incompressible hyperelastic body 
in rotation, by considering a new synthetic and very effective method based 
on complex variables. In doing so, they discovered a striking analogy between 
the equations of motion obtained for a motion general enough to include all 
of the above motions, and the equations obtained for the motion of a nonlin- 



ear string, as considered by Rosenau and Rubin [S3]. The method used by 
Rosenau and Rubin shows in a simple and direct way that all (and more) of 
the different results obtained by Carroll are a direct consequence of material 
isotropy and of the Galilean invariance of the field equations. This is indeed 
an explanation for their ubiquity and versatility. 

The paper is organized in the following manner. In Section II we lay out 
the general field equations of nonlinear elasticity. Because dissipation can- 
not be neglected for most problems in the dynamics of elastomeric materials 
and of soft tissues, we add a simple inelastic tensor of differential type to 
the hyperelastic Cauchy stress tensor. This additional term is similar to the 
classical stress tensor describing a Stokesian fluid; it is the dissipative term 
introduced by Landau and Lifshitz in their monograph on elasticity theory 
P3] and it is commonly used in nonlinear acoustics (see Norris [IZ])- We 
point out that this term is different from the one used in continuum mechan- 
ics as a first approximation of dissipative effects or as a regularization of the 
hyperbolic equations of nonlinear elastodynamics for example in numerical 
computations [H]. Then we specialize the equations of motion to the case 
of one longitudinal wave and two transverse waves. In Section III we seek 
solutions in time/space separable form and derive nonlinear systems of or- 
dinary differential equations for generalized oscillatory shearing motions and 
for generalized sinusoidal standing waves. In Section IV, different classes of 
materials are considered (Hadamard, Blatz-Ko, power-law, fourth-order elas- 
ticity) and some exact solutions are provided. Concluding remarks are made 
in Section V, which recaps the main results. In particular, it is emphasized 
there that a weak nonlinear elasticity theory (up to the fourth order) cannot 
account for some wrinkling phenomenon found in the fully nonlinear theory 
(power- law strain energy). 

2 Preliminaries 

2.1 Equations of motion 

Consider a hyperelastic body with strain energy density S. Let the initial 
and current coordinates of a point of the body, referred to the same fixed 
rectangular Cartesian system of axes, be denoted by X and x, respectively. 
Hence a motion of the body is defined by 

x = x(X,t). (1) 



The response of a homogeneous compressible isotropic elastic solid to de- 
formations from an undistorted reference configuration is described by the 
constitutive relation [3], 

where T-^ is the (elastic) Cauchy stress tensor, 1 is the unit tensor, B is the 
left Cauchy- Green strain tensor defined by 

B := FF^, (3) 

F := (9x/(9X being the deformation gradient tensor, and Ji, I2, I3 are the 
first three invariants of B, 

Ji:=trB, l2:=l[l!-tT{B% J3:=detB. (4) 

To describe the simultaneous effects of thermal and viscoelastic dissipa- 
tion, we introduce the following viscous-like stress tensor j33], 

T^ = 2r/[E - |(tr E)l] + (C + x)(tr E)l, (5) 

where E is the Green-Lagrange strain tensor, 

E:=1(F^F-1), (6) 

and the dot indicates the time derivative. In (|5]), the constant 77 (> 0) is 
the shear viscosity coefficient, the constant ( (> 0) is the bulk viscosity 
coefficient, and the constant x (> 0) is the coefficient related to the ther- 
mal properties of the solid (that is, ambient temperature, thermal expansion 
coefficient, and specific heat). 

Let p and po denote the mass densities of the body measured in the 
current configuration and in the reference configuration, respectively. Then 
the equation of motions, in the absence of body forces, are: div (T^-|-T^) = 
px in current form, or equivalently, 

Div [v^(T^ + T^)(F-^)^] = Pox, (7) 
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in referential form. 



Remark on the dissipative stress tensor 



For incompressible solids, det F = 1 at all times (so that /a = 1 at all times) 
and an arbitrary spherical pressure term, to be determined from initial and 
boundary conditions, is introduced so that 



T^ + T 



D 



as 

-pi + 2— B 



d^^ 1 
oh 



2r7E, 



(8) 



where p = p(x, t) is a Lagrange multiplier. We note that in order to study 
the propagation of finite amplitude motions in viscoelastic incompressible 
materials, some authors [HI [351 EH ES] chose to add to the elastic Cauchy 
stress tensor a viscous term linear in the stretching tensor D, which is the 
symmetric part of the velocity gradient tensor FF~^. In other words, they 
chose to write 



T^ + T 



D 



oh 



2-^B-i + 2z/D, 
dh 



(9) 



for some Lagrange multiplier p^ and for some constant z/. However it must 
be recalled (e.g. Chadwick [2ll p. 146]) that E = F^DF 7^ D in general. 
Thus we remark, for the time being, that the assumptions ^ and ^ are 
not equivalent a priori. We discuss these issues further at the end of the next 
Subsection. 



2.2 Finite amplitude longitudinal and transverse plane 
waves 



Now we consider the following class of motions, 

x = u{X,t), y = Y + v{X,t), z = Z + w{X,t), 



(10) 



which describes the superposition of a transverse wave, polarized in the (YZ) 
plane and propagating in the X-direction with a longitudinal wave propa- 
gating in the X-direction. Then, 
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Here and hereafter, a subscript letter denotes partial differentiation (thus 
ux = du/dX^ Vtt = d'^v/dt'^, etc.) It follows from the definitions ^ that the 
first three invariants of B are given by 



h 



2 + uj^ 



+ 4 



w 



X) 



1 + 2m^ 



+ ^X+Wx, 



and consequently that the strain-energy density S 
function of My and v% + w\ alone, 



^{Ux.V^+W 



x> 



h = u'x, (13) 
S(/i,/2,/3) is here a 

(14) 



With the expressions (11) and (12) substituted into the Cauchy stress ^ 
and into the viscous-like tensor (|5]), the equation of motions ^ reduce to 



Poutt = [{Qi + Q2)ux + \{C + x+ t'n)iu\ + v\ + W\)t]x, 
PoVtt = [Qivx + VVxt]x, 
PoWtt = [QiWx + Wxtlx, 



(15) 



where the functions Qi 
defined by 



Qi{uj^,vl + wl] 



and Q2 



Q2{u\,v\ + w\) are 



Qi 






Qi 






(16) 



Eqs. (15) form a system of two coupled nonlinear partial differential equa- 



tions, generalizing the system derived by Carroll in [9j for an elastic com- 
pressible material. 

The mathematical treatment of initial-boundary- value problems for equa- 



tions such as (15) is not trivial at all. There is a considerable literature on 



the subject and the most updated reference is the recent paper by Antman 
and Seidman [2], to which we refer for further information. Here we use a 
semi-inverse method and we shall consider the possibility to use our explicit 
results to solve initial-boundary- value problems only a posteriori. Obviously, 
the possibility of general uniqueness theorems (such as those established in 
[2]) is quite valuable to give a clear and definitive status of the solutions 
found by a semi-inverse method. 



Now we differentiate (15) with respect to X, and we recast the resulting 
equations in the form 

PoUtt = [{Qi + Q2)U]xx 

PoVtt = [QiV]xx + V^xxt, 

PoWtt = [QiW]xx + riWxxu (17) 

where we have introduced the new functions: U := Ux, V := Vx, W := Wx- 
Further, these equations can be rewritten in an even more compact form, 
by introducing the complex function A, with modulus fl and argument 6, 
defined by 

A(X, t) = n{X, t)e*^(^'*) :=¥ + iW, (18) 



so that 



V = ?ft{A)=Qcos9, W = '^{A) = nsm9, 

V^ + W^ = Q\ (19) 



Then the three equations (17) are 



PoUtt = [{Qi + Q2)U]xx + UC + X+ IvW + ^']xxt, 

PoA« = [QiMxx + V^xxt- (20) 



Now also, the invariants Ji, I2, I3 found in (13) are written as 

h=2 + U^ + n^, /2 = 1 + 2U^ + Q^, h = U\ (21) 

so that Qi = Qi{U'^,VL^) and Q2 = Q2{U'^,^'^) or, more explicitly. 

These equations have been developed in the purely elastic case of a single 
shear wave by Carman and Cramer pj , who determined some numerical and 
asymptotic solutions. 



Remark on finite amplitude motions in incompressible materials 



In incompressible solids, det F 
Then we find 
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It then follows that 

Div [{-pi + 27]E){F^^f] = 
where p := p — {vx + Wx)t, and that 



-Grad p + rjvxxtj + rjWxxtK 



Div [(-pa + 2z/D)(F 



-l\Ti 



-Grad p^ + uvxxtj + uwxxt^- 



(24) 



(25) 



(26) 



Now we see that the equations of motion ([T]) are the same whether the choice 
dsl) or the choice (tol) is made. The two approaches in modeling the dissi- 
pative effects are reconciled, at least as far as finite amplitude motions in 
incompressible materials are concerned. 



3 Solutions in separable form 



The general solution of ( |20[ ) may be found only via numerical methods. To 
obtain some simple analytical information here, we look for reductions of 



the governing equations (20) to a set of ordinary differential equations. The 



underlying idea is to search for A(X, t), defined in (18), in a separable form 
such as 

fi(x,t) = fii(x)fi2(t), eix,t) = ei{x) + e2{t), (27) 

where Qi, 9i are functions of space only and ^2; ^2 are functions of time only. 
Some lengthy computations, not reproduced here, show that in only two 



general cases can equations (20) be separated in such way for any material 
response. 



The results of this Section generalize the results of GarroU [TOl [13] by 
considering more general types of displacement fields and by including vis- 
coelastic effects. We also acknowledge that the solutions presented here are a 
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generalization of the non-linear three-dimensional motions of an elastic string 
found by Rosenau and Rubin [53] (see [25] for a different generalization, to 
finite amplitude waves in rotating incompressible elastic bodies.) 

3.1 Generalized oscillatory shearing motions 

The first general type of separable solution is of the form, 

A(X, t) = [^{t) + i0(i)];te'('=^-^W), (28) 

where fc is a constant and ip, 0, 6 are arbitrary real functions of time. Here 



A is indeed of the separable form given by (18), with 

fii(X) = k = const., n2{t) = [(j){tf + tpity]^, 

ei{X) = kX, d2{t) = d{t) + tan-^[0(t)/V^(t)]. (29) 

Separating real and imaginary parts, we find that the transverse displacement 



field (f,w) derived from (28) is 

v{X, t) = (f){t) cos{kX - e{t)) + ij{t) sm{kX - 9{t)), 

w{X, t) = (j){t) sm{kX - e{t)) - i){t) cos{kX + e{t)). (30) 

In this case, 

n^ = e[(j)''{t) + ij^it)], (31) 

is a function of time only. 

Now the first equation of motion (pOJ)i reduces to 



PoUtt = [{Qi + Q2)U]xx + UC + X + |^)[f/V«, (32) 

Qi = Qi(f/2(X,t),fi2(t)), Q, = Q,(U\X,t),n\t)). (33) 



On the other hand, the second equation of motion (20)2 reduces to 



Po Ait = Q 1 Axx + V^xxt ■ (34) 



Using the separable form (28) for A, factoring out the exponential term, and 



separating the real and imaginary parts yields 

Po{^9" - i^e'^ + 2^'9' + ^") + eQi^ + er]{4>9' + ^') = 0, 

Po(0" - 24j'e' - (pe'^ - tpe") + eQ^<i) + k^r]{(p' - ^e') = o. (35) 
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This is a nonlinear system of two ordinary differential equations in the three 
unknowns (f){t), 0{t), and tp{t) either: when Qi is a constant independent 

= Qi(t). By (Qi, this 
U{t), a function of time 



of its arguments f/^(X, t) and f2^(t), or: when Qi 



latter condition is satisfied if and only if U{X, t) 

only. It then follows by (33)2 that Q2 = Q2{t) also. By substitution into the 

first equation of motion (20 )i, we obtain PqU" = 0. Then the longitudinal 



displacement field is given by 



U{t) = C^t + C2 = ux, 



(36) 



so that 

uiX, t) = {Cit + C2)X + C3, (37) 

where Ci, C2, C3 are constants. This is a homogeneous accelerationless mo- 
tion. 

Hence we showed that for any compressible elastic material with viscoelas- 
tic dissipative part defined as in ([5]), maintained in a state of longitudinal 
simple extension, the special transverse waves (30) may always propagate. 



Because (35) is a system of only two equations in three unknowns, there 



is a certain freedom in considering special classes of solutions. One such class 
of special solutions is obtained by considering ?/^ := in (30), leading to 



v{X,t) =(l){t)cos{kX-9{t)), 
w{X,t) =(j){t)sm{kX-e{t)), 



(38) 



a direct generalization of the classical damped harmonic circularly-polarized 
wave solution, 



v{X, t) = Ae-''* cos{kX - ut), 
w{X, t) = Ae-''^ sm{kX - uot), 



(39) 



where A, h, and u are suitable constants. For these special motions at ■i/' := 0, 



the system (35) reduces to 



Po{(j)e" + 2(j)'e')+r]k\ 
po(0" - 0^") + k^Qi4 



rjk 



0, 

2 J,/ 



0. 



The equation (40)i admits the first integral 

6^9' = Ee "0 



(40) 



(41) 
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where E is a constant of integration. 

If ii^ = 0, then 6{t) = Oq, a constant. We end up with standing waves, 

v{X,t)=(j){t)cos{kX -Oo), w{X,t) = (f){t)sm{kX -Oo), (42) 



where, according to (40)2, is a solution to 

+ k^7](j)' + eQi {k^(j)^, t)0 = 0, (43) 



which is the equation of a damped nonhnear oscillator. 



If i? 7^ 0, then solving (41) for 9', we end up with the single equation 



(40)2, which reduces to 



+ A;V + k'Qi(j) - ^e"^* = 0. (44) 

In the elastic case [r] = 0), this equation describes the plane motion of a 
particle in a central force field. In general [r] ^ 0), it is a nonlinear and 
non-autonomous differential equation. 



For harmonic wave propagation, 6{t) is taken as 6{t) = cut. By (41) 



is proportional to e ^po and describes the usual linear damping function. 



However by substitution in (40)2, we see that this situation is possible only 



if the material response is such that Qi is a constant. 



In this Subsection we analyzed the solutions (38) which are the analogue 
of the Class I solutions presented by Rosenau and Rubin [31] for purely elastic 
strings. Qualitatively the solutions describe particles moving in a helical path 
which lies on a cylindrical surface of time- varying radius. 

3.2 Generalized sinusoidal standing waves 

We now consider the second class of solutions in separable form, namely 

A{X, t) = [{^4>{X) + HX))9\X) + (0'(X) - zV^'(X))]e'(-*+^(^)), (46) 
where A; is a constant and ipiX), (f){X), 0{X) are arbitrary functions of space. 



Here A(X, t) is indeed of the separable form (18), with 

ni{X) = [{(j)' + 4)9' f + {(1)9' - ^')li n2{t) = 1 = const., 

9i{X) = 9 + tan-i[(0^' - ij')/{^9' + 0')], ^2(t) = ujt. (47) 
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The transverse displacement field {v,w) corresponding to (46) is 



v{X, t) = 0(X) cos(wt + 9{X)) + iP{X) sm{ujt + e{X)), 

w{X, t) = 0(X) sin(wt + e{X)) - ^(X) cos{iot + e{X)). (48) 

In this case, Q = Qi{X), a function of X only, and the first equation of 

(49) 



motion (20)i reduces to 

PoUtt = [(Qi + Q2)U]xx + U( + x + lv)[U^]xxt 
where 



(50) 



Now, introducing (46) into the second equation of motion (20)2 and sep 



arating the real part from the imaginary part, we obtain 

- p^uj^i^ff + 00 = {i^ff + <P'){Qi)xx + 2(2^'^' + ^r + r - <pe'^){Qi)x 

- r]uj{?><p"e' + 30'^" + 0^'" - (pe'^ + 3^'e'2 + ?>0e" - ^"'), (si) 

and 

- p^u\<pe' - ^') = {<pe' - i^'){Qi)xx + 2(20'^' + (t>e" - v^" + ¥^){Qi)x 
+ (30"^' + 30'^" + (pe'" - 0^'^ + ?,^ij'e"^ + 3^0'^" - ^"')Qi 
+ T]uj{?>i)"e' + 3V^'^" + i^e'" - ^^'3 - 30'^'' - ?,<pe'e" + 0'"). (52) 



A sufficient condition for the equations (51) and (52) to compose a nonlinear 



system of two ordinary differential equations in the three unknowns (p , ip, 
and 6 is that the longitudinal field depends only on X: U{X,t) = U{X). In 



this case, Qi = Qi{X), Q2 = Q2{X) and (49) reduces to 

[{Qi{X) + Q2{X))U{X)]" = 0. (53) 

Further progress is made by fixing one of the three unknown functions. 



9 := 0, say. Then the governing equations (51) and (52) reduce to 

(<5i0')' + Pow^0 = -r]^i^", iQii'')' + Po^^i' = r]^4>"- (54) 

At ri = 0, these equations are formally equivalent to the equations derived 
by Carroll [H], with the difference that the materials response functions Qi 
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and Q2 depend not only on Vl? (the amount of shear) but also on U . In 
Carroll [15] (see system 5.12 in that reference), the attention is restricted to 
incompressible fluids and solids, for which U{X) must be constant. 

The solutions considered in this sub-Section contain the Class II solutions 
of Rosenau and Rubin 



4 Specific materials 

Now that we have established the possibility of reducing the general nonlin- 
ear partial differential equations to ordinary differential equations, we must 
specify constitutive relations in order to make progress. 

4.1 Hadamard materials 

Finite amplitude motions in elastic Hadamard materials have been thor- 
oughly studied, see John |12] and Boulanger et al. [fij. The analysis con- 
ducted in this paper allows for a treatment of a dissipative Hadamard solid. 
The Hadamard strain energy function is defined by 

2S = C{h - 3) + D{h - 3) + G{h), (55) 

where G(/3) is a material function and C, D are material constants such 
that [6]C>0,D>0, orC>0,D>0. The connection with the Lame 
constants A and fi of the linear theory of isotropic elasticity is made through 
the relations: 

C = 2/i + G'(l), D = -fi-G'{l), 2G"{l) = X + 2fi. (56) 



For this material, the functions Qi and Q2 defined in (16) are 

Qi = C + D, Q2 = D + G'{U^), (57) 



and the equations of motion (20) reduce to 



PoUu = {[G + 2D + G'iU^)]U}xx 

+ UC + x + lv)iu' + n')txx, 

PoAtt = {G + D)Axx + r]Atxx. (58) 



14 



We point out that the equation governing the propagation of the trans- 



verse waves (58)2 is always hnear, for all Hadamard materials, elastic or 



dissipative. Hence we see at once that classical harmonic (damped and at- 
tenuated) transverse waves are always possible for these materials. 

The material function ^(/s) accounts for the effects of compressibility. 
For example, Levinson and Burgess |11] proposed the following explicit form. 



Gih) = (A + fi){Is - 1) - 2(A + 2/i)( V/3 - 1). (59) 

This function leads here to a remarkably simple system of equations: 

PoUtt = (A + 2fi)Uxx + U( + X+ lv){U^ + fi')txx, 

PqAu = pAxx + V^txx- (60) 

Now we take a look at some of the solutions investigated in the previous 



Section. Classical damped harmonic circularly polarized waves such as (39) 



are possible in a dissipative Hadamard material because then fl = kAe , a 
function of time only, and Qi = C + D = fj,, a. constant. The corresponding 



dispersion equation (45) is here 



00 
P = Po— + 






4/00 



(61) 



The associated longitudinal motion may be either the homogeneous acceler- 



ationless motion (|36|) or any solution to the nonlinear equation (|58|)i, here: 

(62) 



PoUu = {[G'{U') - G\l)]U}xx + U( + X + lv){U')txx, 
which for the Levinson and Burgess choice ( [59| of G is: 

PoUu = (A + 2fi)Uxx + U( + X + IvW^xx. 



(63) 



Now consider the sinusoidal standing waves of Eq.([48|) at 0{X) = 0. 
When the longitudinal field U depends on X only, their behavior is governed 



by Eqs.(54), here: 



/i0" = pouP'ct) + rjOJip" , pip" = PoOj'^'ip — V'^'P" J 
or equivalently, by the single complex equation 

(/i — iriuj){(j) + ii^)" = P^'^i'P + 'i'i')- 



(64) 



(65) 



15 



This equation possesses the following class of attenuated solutions, 

-"^[h sin(/3X) + k2 cos(/3X)], 



0(X) 
^|J{X) 



e""^ [A;3 cos(/3X) + k^ sin(/3X)], 



where ki, k2 are disposable constants and 

a, (3 = 



pouj^ 



2(^2 + ^2^2 



V /i^ + r/2|^2 _|_ ^ 



(66) 



(67) 



These solutions may be used to solve some simple boundary value problems. 
For example, consider the case of semi-infinite body bounded by a vibrating 
rigid plate at X = 0, and take the velocity components of the plate as 

a; = 0, y = Vucos{ut), z = Vcosm^ut). (68) 

Then 0(0) = 0, 4'{0) = V, and the associated transverse displacements are 

-°^[sin(/3X) cos(a;t) + Vcos{/3X) sin(a;t)]. 



viX.t)-- 
w{X,t) 



-aX 



sm{f3X) sm{ut) - V cos(/3X) cos(a;t)]. 



(69) 



The associated longitudinal displacement independently satisfies (53), here: 

{[G'{U^)-G'{1)]U}" = 0, (70) 



which for the Levinson and Burgess choice (|59J) of G is simply: U" = 



Then the basic displacement (69) may be superimposed upon an axial static 



stretch or a homogeneous motion in the axial direction. In a similar way, 
it is possible to solve the case of a slab fixed at X = and oscillating at 
X = L, or it is possible to use the various integration constants to fix the 
stress traction or shear stress at the boundary of a vibrating half-space or of 
a slab. 

4.2 Blatz-Ko materials 

Experiments on compressible polyurethane rubber lead Blatz and Ko [3] 
to propose some strain energy functions which have since received much 
attention, see for instance, Beatty ^. In particular, two reduced forms of 
the Blatz-Ko general constitutive equation were deemed appropriate to model 
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certain polyurethane rubber samples. One is the Blatz-Ko strain energy 
function for a solid, polyurethane rubber, 



/^ 



1//3 



[h-3 + ^{I^^^ - 1)] 



(71) 



where /x and /3 are constants. Direct comparison with (55) shows that this 
material is a special Hadamard material, with the identifications: C = fi, 
D = 0, G^Is) = —ji(3I^. The other is the Blatz-Ko strain energy function 
for a foamed, polyurethane elastomer, 



2 



^-3 + /3(/V^ 



1) 



(72) 



where /x and /3 are constants. The continuity with linear elasticity requires 
that /i is the infinitesimal shear modulus and that /3 = (1 — 2z/)/z/, where u 
is the infinitesimal Poisson ratio. Blatz and Ko's experiments showed that 
typically, u = 1/4 (and so, (3 = 2). 



For this latter material, the functions Qi and Q2 defined in (|16|) are 



2(/?-l) 



(73) 



and the equations of motion (20) reduce to 



PoUu = -i^ {i + n^)u-''-u 



3/3-2 



XX 



+ UC + x + h){u' + n%xx, 



PoAtt = ij{AU' 



)xx 



+ riA 



txx- 



(74) 



Here, and in contrast with the case of Hadamard materials, the transverse 
wave is coupled to the longitudinal wave. 

First, consider the generalized oscillatory shearing motions of Section 
III.A. As seen there, they are governed by nonlinear ordinary differential 
equations when the longitudinal displacement u is the homogeneous accel- 



erationless motion of (86)2. Then the standing waves of (42) are governed 
by 



" + r]k 



2 J,/ 



0. 



(75) 



At Ci = 0, this is the classical equation of a damped oscillator, whose 
type of motion (exponential and/or sinusoidal) is decided by the sign of the 
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quantity: {rjkC2Y — ^PolJ'] hence the nature of the transverse motion depends 
on the longitudinal extension via the parameter C2. 

At Ci 7^ (linear stretch rate), we perform the change of variables C, = 
Cit + C2, so that (75) reduces to 



^nd^(/) ^^^od^ r 

where the quantities T, V are defined by 

rjk"^ 



0. 



T 



2poCi 



r:= 






(76) 



(77) 



In the purely elastic case (r^ = 0), T vanishes and the solution to this equation 
is quite simple: we find that for 



fciC^+ + ^2C^" , where A± = (1 ± Vl - r)/2; 

A;iv/C + /^2yclnC; 

A;iv^cos(v/r^TlnC) + A;2V^sin(v/r^^lnC); (78) 

where ki and k2 are integration constants. We note that all these solutions 
blow up as t ^ 00 (equivalent to C — ^ 00); that for F < 1, the solutions are 
monotonic increasing; and that for F > 1, the solutions have a slight oscilla- 



<1, 


m 


= 1, 


0(C) 


>1, 


0(C) 



tory character. In the viscoelastic case (r/ 7^ 0), equation (76) is solvable in 



terms of special functions and a richer variety of solutions emerges, because 
now the solutions are not necessarily unbounded. Moreover, although they 
are damped solutions, they do not necessarily vanish with time. In order to 



solve equation ( 76 ) explicitly, we introduce the function / defined by 

0(C) = e-^^/(2TC), 
and we find that it satisfies a Whittaker equation: 



r(2Tc) + 



F 



4 ^ 4(2TC)2 



/(2TC) = 0. 



(79) 



(80) 



It follows that the solution to ( 76 ) is expressed in terms of Whittaker's func- 
tions as: 



0(C) 



-TC 



hW. 



n-r 



(2TC) + k^M^ 



o,|v^(2TC) 



(81) 




Figure 1: Influence of tlie stretcli rate Ci in tlie X-direction on tlie sliearing 
motion solutions for a Blatz-Ko material (foamed, polyurethane elastomer). 

where ki and k2 are integration constants. To illustrate the influence of the 
stretch rate in the X-direction (through the constant Ci) upon the behavior 
of the solutions, we consider that the material is unstretched at t = (so 
that C2 = 1 and 0(t = 0) = 0); we take 0'(t = 0) = 1, fik'^/po = I, V = /^/lO 
(so that T = l/(20Ci) and T = A/Cf); and we let Ci = 0.2, 1, 1.5. Figure 
1 shows that for a "small" rate, the solution undergoes a few oscillations 
before it vanishes completely; as Ci increases, the oscillations disappear but 
the solution tends to an increasing finite asymptotic value. We argue that 
this behavior can be related to the problem of a tensile impact: when the 
material is plucked at "low" speed, then a few oscillations take place before 
the material returns to its original state; when the material is plucked at 
"high" speed, then the material can accomodate an asymptotic transverse 
wave. 

Now, consider the generalized sinusoidal standing waves of Section lll.B. 
As seen there, they are governed by nonlinear ordinary differential equations 
when the longitudinal displacement u is such that ux = U = U{X), a 
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function of X only. This condition leads to (53), here integrated twice as 

(0'2 + ^'2 _ ip-3 + ^^ ^ CsX + C4, (82) 



where C3, C4 are disposable constants. Taking the typical value (3 = 2 
corresponding to Poisson ratio u = 1/4, we find 



U 



(p^ + 4)' 



1 



C3X + C4 - 1 



(83) 



Hence for this material, the governing equations (54) for the transverse so- 
lutions of type (48) at 9{X) = are, using (83) and (73)i, 



/i 



/i 



02 +?/) 



1 



C3X + C4 - 1 

CsX + C4 - 1 



2 -■ ' 
3 



^' 



+ PqU (f) = -r]uip'' 



+ pou'^i' = -qujcj)". 



(84) 



The final governing equations are highly nonlinear, even with the choice 
/3 = 2. The possibility of solving them in closed form seems quite remote. Of 
course, numerical and qualitative analyses may be performed with the usual 
methods of dynamical systems theory. 

4.3 A separable strain energy density 

Recall that the strain energy function of an incompressible material. Sine say, 
corresponds to the restriction to the subspace (Ji, J2, 1) of a strain-energy E 
say, defined in the full space (Ji, /2, h)'- 

Si„c(/i,/2):=S(/i,/2,l). (85) 

Hence the classical incompressible neo-Hookean form, 

Sine = p{h - 3)/2, (86) 

may be associated for example with the full-space strain energy 



E = |(/i-3-21nJ), J:= v^. 



(87) 
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often used in the classical molecular theory of rubber 

Now, constitutive equations for compressible hyperelastic materials come 
in many formulations. One of them consists in adding a purely volumetric 
term Evoi(^) say, to a basic strain energy density function S, whose restriction 



(85) is the strain energy function of an incompressible material T^mc- Then 
the final strain energy density for a compressible hyperelastic material can 
be written as ^ 

S = S(/i,/2,/3) + S™i(J). (88) 

Several models have have been proposed in the literature for this pure 
volumetric part (or bulk term,) of the strain energy function. Ogden [18] 
proposed the form: 

S^„i(J) = A/3-'(/51nJ+J-^-l), (89) 

where A is the first Lame modulus and /3 (> 0) is an empirical parameter. 
Flory [27] proposed the form 

KUJ) = ic/2){\nJ)\ (90) 

and Simo and Pister [56], the form 

KiliJ) = ic/2)[{J'-l)-\nJ], (91) 



which is (89) at /3 = —2. Recently Bischoff et al. |5.j proposed the form 



c 



S;oi(^) = 7^[cosh/3(J-1)-1]. (92) 



(3 

In this manner, several strain energy functions may be constructed. As 
an example, consider the power-law function of Knowles 



C 

J 



l + -(/i-3)) -1 
n 



(93) 



where C, b, and n are constitutive parameters, all three assumed positive. 
This simple strain energy captures several important features of rubber-like 



materials. At n = 1, it recovers the neo-Hookean form (86); at n > 1 it 
describes strain-stiffening materials; at n < 1, strain-softening materials. 
These hyperelastic properties prove crucial to the understanding of complex 
phenomena such as dynamic fracture, as shown recently by Buehler et al [7]. 
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Then, using the formulation exposed above, we may construct the following 
compressible strain energy function 



C 

~b 



l + -(/i-3) 



n 



ClnJ + S™i(J). 



(94) 



Turning back to the finite amplitude motions of Section II, we find that 
J = y/h = ux = U, so that Svoi(^) = Svoi(f^), and the functions Qi, Q2 



defined in (16) are here 



in-l 



Qi = 2C[i + -iu' + n'-i)Y 

Q2 = U-'KoiiU) - CU-\ (95) 

In general, the corresponding equations of motion are coupled and rather 



involved. Consider the special case of the standing waves (42), with the 



choice U{t) = C2 (Ci = in (36)). Then 0(t) satisfies (43), which is here 

0. (96) 



/," + ^0' + 2— [1 + -{cl - 1 + k^^)T-\ 

Po Po n 



We notice that this equation may admit not only the static solution: 
= 0, but also the nontrivial static solutions: 



±il/k)^l-Cl-n/b. 



(97) 



These solutions are real if and only if 6(1 — Cf ) — -n, > 0: when the body is 
compressed in the X-direction {u = C2X, C2 < I), then the solution is real 
if 6 > n/(l - C2) > 0; when the body is stretched {u = C2X, C2 > 1), the 
solution is not real (recall that Knowles [10] assumed that 6 > in (93)). 



Another important general property of equation (96) is that it can be 
derived [5^ from Lagrangian 



fc^i7 



/:=le^*[0'2_v(0)] 



(98) 



where 



V(0) 



ec 



po 



[1 + -{C; 



1 + A:^C')]""'CdC 



C 
Pob 



[i + -{ci-i + k'<j)')r. 



n 



(99) 
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As an example we consider a strain- stiffening material, by taking n = 2 



in (94). Then (96) is a Duffing equation with damping, 

0" + 6(f)' - /30 + a(f)^ = 0, (100) 

where 

6:=^, /3:=^^[6(l-C2)-2], a:=b^^. (101) 

Po Po Po 



Clearly, a > 0, 6 > 0. To make the connection with results by Holmes [37], we 
impose /3 > also, which happens only when the body is compressed in the 
X-direction {u = C^X, Ca < 1) and when h > 2/(1 - C|). As seen for ( [TOOJ ), 



this equation possesses three fixed points in the phase plane, namely: (0,0), 
which is a saddle, and {±^/pJa, 0), which are two sinks. Holmes |37j showed 



that equation (101) is locally and globally stable. It follows that as t — )■ oo. 



our standing waves must approach one of the four wrinkled configurations 
\imv{X,t) = ±^/l3jacos{kX -60), 
\imw{X,t) = ±^/JJ^sm{kX-9o), (102) 



t—^00 



for almost any initial conditions. As any one of the four final configurations 
is equivalent to another, we conclude that the material is destabilized by 
the waves (only some special choices of the initial conditions will lead to the 
unstressed deformation corresponding to the saddle point (0,0); these special 
initial conditions are found by computing the corresponding stable manifold 
in the phase plane, and we refer to Holmes |37] for the details.) 

The situation described here can be extended to all positive integers n, 



because the leading of nonlinearity of the potential V(0) in (99) is 0^" so 
that the leading order of nonlinearity in the resulting differential equation 
is always to an odd integer power. We also point out that for a given b, an 
increasing n requires a decreasing C2 for the appearance of the fixed points 



(97) that is, a greater compressive stretch in the X-direction. This is in 
agreement with the "physical" expectation that the stiffer a material is, the 
harder it is to destabilize it. 

For non-integer powers n, the situation is more complex because destabi- 
lizing configurations may or may not appear in the case of strain-hardening 
materials (n > 1) and in the case of strain-softening materials {n < 1). 
Obviously the dynamical systems uncovered here are susceptible to a more 
complete and more careful analysis than the one provided, which may reveal 
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even richer behaviors than those already described, such as for example the 
appearance of parametric resonance; these aspects are however beyond the 
scope of the present paper. 

Finally, we note that the instability evoked in this Subsection is somewhat 
reminiscent of the buckling of an elasticita under a compressive load, but we 
point out several differences. The buckling of an elastica is a phenomenon 
related to geometrical non-linearities, whereas our instability is due to con- 
stitutive non-linearities; in an elastica, more than one inflexion may occur 
whereas here, we could multiply instabilities only by considering a multiple 
well strain energy function (and such a strain energy is not usual for the 
modehng of soft tissues and elastomers); for an elastica, the stabihty of the 
various fixed points must be examined thoroughly, whereas here, the stability 
of such points comes out directly from our dynamical approach. 

4.4 Fourth-order elasticity 

In the linearized theory of "small-but-finite" amplitude waves, the strain 
energy function must be expanded up to the fourth order in the strain, in 
order to reveal nonlinear shear waves ^9]. In that framework, Murnaghan's 
expansion |16] is often used (see for instance Porubov [52]): 

^ A + 2/i .2 ^ . / + 2m o ^ . . . .4 .n. . . .2 

L = 1^ - 2fll2 H li - 2miil2 + ms + Vil-^ + Z/2Zi«2 + ^^S^l^S + ^^4^2; 

(103) 
where A, /x are the Lame moduli, /, m, n are the third-order moduli, and 
T^ii ^'2, t'a, 1^4 are the fourth-order moduli (other notations exist: see Norris 
[17j for the connections). In this expansion, we used the first three principal 
invariants zi, 12-, is of E, the Green-Lagrange strain tensor; they are related 
to the first three principal invariants Ji, I2, I3 of B defined in Q by the 
relations [1], 

h = 2h + 3, 

12 = Ah + 4^2 + 3, 

13 = 2ti + 4^2 + 8z3 + 1. (104) 

By chain rule differentiation, we find that the functions Qi and Q2 defined 
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in (16) are here, 



Q^ = ;i + (A + 2/i + m)ii + (/ + 2m - 



^^2^2 



(2m + z/4)i2 + 4:iyiil + 2z/2iii2 + ^^3^3, 



Z/2 2 
Q2 = -fJ'- mil + — 2^ + Z/4«2, 



where now zi, Z2, ^3 are found from (|104) and (21) as 



«i = ([/' + fi' - l)/2, 



«2 



-fiV4, ^3 = 0. 



(105) 



(106) 



Several comments are in order at this point. First of all, we recall that 



the equations corresponding to ( [15| ) for third-order elasticity (without dissi- 

in 1960. 



pation) were first derived by Gold'berg 

Then we note that recently there has been a renewed interested in fourth- 
order elasticity following experiments on the nonlinear acoustic properties of 
soft tissue-like solids [20]. Hence Hamilton et al. [32j proposed the following 



reduced version of the expansion (103), suitable when the relative portion of 



energy stored in compression is negligible, as is expected for shear deforma- 
tions and motions of soft tissues, 



-2/iZ2 + m3 + 1/4^2. 



In this case, the functions Qi and Q2 defined in (16) simplify to 

Ql = -Q2 = At - ^^4^2 = /i + I^4fi^/4, 



(107) 



(lOJ 



and we check at once that the longitudinal equation of motion (15)i becomes 
a trivial identity in the purely elastic case (C = X = ''7 = 0), because here 
Ql + Q2 = 0. 

Finally, we also notice that the theories of third- and fourth-order elas- 
ticity have been used to study bulk solitons propagation in elastic materials. 
For instance, Hao and Marris [33] discuss the issue of acoustic solitons from 
both experimental and theoretical points of view; they consider the possibil- 
ity to produce KdV solitons by adding a fourth-order spatial derivative to the 
longitudinal wave equation and by performing some ad hoc approximations. 



The general expressions (103), (105), and (106) yield quite complicated 



equations of motion. If for example we restrict our attention to the standing 



waves (42), we find that the polynomial nonlinearity of the resulting ordinary 



25 



differential equation is of the fifth order for third-order elasticity, and of 
seventh order for fourth-order elasticity. If we further restrict our attention to 



the reduced expansion (107) we find that equation (43) reduces to a damped 
Duffing equation, 

0" + ^0' + ^(4/i + ^2^,^02)0 = 0. (109) 

This equation is characterized by a positive linear stiffness and therefore 
the peculiar behavior found in the previous sub-section for some power-law 
materials is ruled out. For completeness, we point out that Holmes and Rand 



computed an approximate solution of equation (109), using the method 
of averaging. 

5 Concluding Remarks 

We studied the propagation of finite amplitude waves in a general nonlinear 
elastic material, with dissipation taken into account by means of a simple 
mechanism of differential type. We showed that the general equations of mo- 
tion admit some beautiful and most interesting reductions to ordinary differ- 
ential equations, using a direct method based on complex functions. Clearly, 
these reductions are a consequence of the inherent symmetrical structure of 
the balance and constitutive equations and may therefore be recovered also 
by the standard methods of group analysis; we argue however that our direct 
method is more simple and more revealing from a mechanical point of view. 
Moreover, our results complement and generalize the well-known results of 
Carroll in various directions, as already pointed out in the Introduction; 
they also emphasize the analogy between Carroll's solutions and some exact 
solutions proposed for nonlinear strings by Rosenau and Rubin [[53] . 

The solutions provided here can be used not only as benchmarks for the 
more complicated numerical analysis required in "real life" applications, but 
also for a better understanding of the mechanical properties and of the mathe- 
matical structure of various usual models. These advantages were highlighted 
in the discussion on standing waves solutions conducted in Section III. There, 
we focused much of the discussion on the special case of standing waves for 
the sake of simplicity and brevity, but it is clear that our general methods 
of investigation can be applied to all the ordinary differential equations that 
were derived. 
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One of the important findings uncovered by the analysis of our solutions 
is the possibility to predict the appearance of highly symmetric wrinkles in 
an elastic medium compressed longitudinally. This phenomenon has been 
observed experimentally by several researchers when, for example, a hard 
film is deposited on a soft material and put under compressive strain [5U] . 
Our analysis showed that these wrinkles may indeed appear when the com- 
pression parameter satisfies some simple inequalities; it also showed that the 
detailed viscoelastic behavior of the material is unimportant because of the 
highly attracting nature of the wrinkled states of deformations. We note of 
course that our solutions have a higher symmetry than the patterns usually 
observed in experiments, probably because our solutions are characterized 
by strong invariance and because they are bulk solutions, in contrast to the 
half-space geometry generally used in experiments. Despite these limitations, 
we found that these solutions are compatible with the relatively simple non- 
linear power-law strain energy density whereas weak non-linear theories (up 
to the fourth-order) are not adequate to predict these features. 

Our approach allowed a direct and complete comparison of the models 
used in nonlinear acoustics and in continuum mechanics. It was seen that 
although the popular weak nonlinear theories are often used to study some 
special features of wave propagation, a consequence of their polynomial na- 
ture is that they cannot contain all the necessary information of the general 
constitutive theory. 

Last but not least, we gave some explicit solutions to some models of 
compressible hyperelasticity. New exact solutions are always welcomed and 
valuable additions to the short list of exact solutions found in the literature. 

We conclude by pointing out that several generalizations of our results 
are possible in principle. For instance, the consideration in polar coordinates 
of a wave propagating in the radial direction coupled to two shear waves in 
the axial and azimuthal direction is one of such possible generalizations. In 
this case the following class of motions is considered 

r = r{R,t), e = e + g{R,t), z = Z + u{R,t), (110) 

where {R, G, Z) are cylindrical polar coordinates associated with the unde- 
formed configuration and (r, 6, z) are cylindrical polar coordinates associated 



with the deformed configuration. It is possible to show that from (110) and 
the corresponding kinematical quantities of interest, the equations of motion 
may be reduced to a set of three partial differential equations, similar to 
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what we did for the motions (10) in (15). Unfortunately, it seems that in 



this case the equations of motion do not allow general solutions in separable 



form, as was the case in Section III for the equations (15)). Of course, we do 
not exclude that for special classes of materials, it is possible to find some 
classes of wave solutions (see for example [21] for compressible materials and 
for incompressible materials). 



References 

[1] S.S. Antman, Nonlinear Problems of Elasticity (Springer, New York, 
2004). 

[2] S.S. Antman and T.I. Seidman, The parabolic-hyperbolic system gov- 
erning the spatial motion of nonlinearly viscoelastic rods. Arch. Ration. 
Mech. Analysis 175 (2005) 85-150. 

[3] P.J. Blatz and W.L. Ko, Application of finite elasticity to the deforma- 
tion of rubbery materials. Trans. Soc. Rheol. 6 (1962) 223-251. 

[4] M.F. Beatty, Topics in finite elasticity: Hyperelasticity of rubber, elas- 
tomers, and biological tissues - with examples, Appl. Mech. Rev. 40 
(1987) 1699-1733. 

[5] J.E. Bischoff, E.M. Arruda, and K. Grosh, A new constitutive model for 
the compressibility of elastomers at finite deformations. Rubber Chem. 
Tech. 74 (2001) 541-559. 

[6] Ph. Boulanger, M. Hayes, and C Trimarco, Finite-amplitude waves in 
deformed Hadamard elastic materials, Geophys. J. Int. 118 (1994) 447- 
458. 

[7] M.J. Buehler, F.F. Abraham, and H. Gao, Hyperelasticity governs dy- 
namic fracture at a critical length scale. Nature 426 (2003) 141-146. 

[8] G.P. Carman and M.S. Cramer, Nonlinearity parameters for pulse prop- 
agation in isotropic elastomers, J. Acoust. Soc. Am. 91 (1992) 39-51. 

[9] M.M. Carroll, Some results on finite amphtude elastic waves. Acta Mech. 
3 (1967) 167-181 . 



28 



[10] M.M. Carroll, Oscillatory shearing of nonlinearly elastic solids, Z. angew. 
Math. Phys. 25 (1974) 83-88 . 

[11] M.M. Carroll, Plane elastic standing waves of finite amplitude, J. Blast. 
7 (1977) 411-424. 

[12] M.M. Carroll, Finite amplitude standing waves in compressible elastic 
solids, J. Blast. 8 (1978) 323-328 . 

[13] M.M. Carroll, Reflection and transmission of circularly polarized elastic 
waves of finite amplitude, J. Appl. Mech. 46 (1979) 867-872. 

[14] M.M. Carroll, Unsteady liomothermal motions of fluids and isotropic 
solids. Arch. Ration. Mech. Analysis 53 (1974) 218-228. 

[15] M.M. Carroll, Plane circular shearing of incompressible fluids and solids, 
Q. J. Mech. Appl. Math. 30 (1977) 223-234. 

[16] J.H. Curran and M.M. Carroll, Sinusoidal shearing of liquid crystals. 
Int. J. Bng. Sci. 18 (1980) 971-977. 

[17] M.M. Carroll, On circularly-polarized nonlinear electromagnetic waves, 
Q. Appl. Math. 25 (1967) 319-323. 

[18] M.M. Carroll and M.F. McCarthy, Finite amplitude wave propagation in 
magnetized perfectly electrically conducting elastic materials. In: Bias- 
tic Wave Propagation., M.F. McCarthy and M.A. Hayes, eds. (Elsevier, 
North-Holland, 1989) pp. 615-621. 

[19] C. Carstensen and G. Dolzmann, Time-space discretization of the non- 
linear hyperbolic system Uu =dw{(T (Du) + Dut) , SI AM J. Numer. 
Analysis 42 (2004) 75-89. 

[20] S. Catheline, J.L. Gennisson, and M. Fink, Measurement of elastic non- 
linearity of soft solid with transient elastography, J. Acoust. Soc. Am. 
114 (2003) 3087-3091. 

[21] P. Chadwick, Continuum Mechanics (Dover, New York, 1999). 

[22] P.J. Chen, Growth and decay of waves in solids, in: C.A. Truesdell (ed.) 
Handhuch der Physik, Vol. VIa/3 (Springer, New York, 1973). 



29 



[23] CM. Dafermos, Hyperbolic Conservation Laws in Continuum Physics, 
(Springer, New York, 2000). 

[24] M. Destrade and G. Saccomandi, Finite amplitude inhomogeneous waves 
in Mooney-Rivlin viscoelastic solids. Wave Motion 40 (2004) 251-262. 

[25] M. Destrade and G. Saccomandi, Some results on finite amplitude elastic 
waves propagating in rotating media. Acta Mech. 173 (2004) 19-31. 

[26] B. Erman, Molecular aspects of rubber elasticity, in: G. Saccomandi 
and R.W. Ogden (eds.) Mechanics and Thermomechanics of Rubberlike 
Solids, CISM Lecture Notes 452 (Springer, New York, 2004). 

[27] P.J. Flory, Thermodynamic relations for high elastic materials. Trans. 
Faraday Soc. 57 (1961) 829-838. 

[28] R.L. Fosdick and G.P. Macsithigh, Shearing motions and the formation 
of shocks in an elastic circular tube, Q. Appl. Math. 38 (1980) 191-207. 

[29] Z.A. Gold'berg, Interaction of plane longitudinal and transverse wave, 
Akusticheskii Zhurnal 6 (1960) 307-310. 

[30] J. Hadamard, Legons sur la Propagation des Ondes et les Equations de 
I'Hydrodynamique (Hermann, Paris, 1903). 

[31] J.B. Haddow and R.J. Tait, A note on the propagation of finite am- 
plitude shear waves on a hyperelastic solid, ASME J. Appl. Mech. 53 
(1986) 470-471. 

[32] M.F. Hamilton, Y.A. Ilinski, and A. A. Zabolotskaya, Separation of 
compressibility and shear deformation in the elastic energy density, J. 
Acoust. Soc. Am. 116 (2004) 41-44. 

[33] H.Y. Hao and H..J. Marris, Experiments with acoustic solitons in crys- 
talline solids, Phys. Rev. B 64 (2001) 064302. 

[34] M.A. Hayes and G. Saccomandi, Finite amplitude transverse waves in 
special incompressible viscoelastic solids, J. Blast. 59 (2000) 213-225. 

[35] M.A. Hayes and G. Saccomandi, Finite amplitude waves superimposed 
on pseudoplanar motions for Mooney-Rivlin viscoelastic solids. Int. J. 
Nonlinear Mech. 37 (2002) 1139-1146. 

30 



[36] M.A. Hayes and G. Saccomandi, Antiplane shear motions for viscoelastic 
Mooney-Rivlin materials, Q. J. Mech. Appl. Math. 57 (2004) 379-392. 

[37] P. Holmes, A nonlinear oscillator with a strange attractor, Phil. Trans. 
Roy. Soc. London A 292 (1979) 419-448. 

[38] P. Holmes and R. Rand, The bifurcations of Buffing's equation: an 
application of catastrophe theory, J. Sound Vibr. 44 (1976) 237-253. 

[39] Z. Huang, W. Hong, and Z. Suo, Evolution of wrinkles in hard films on 
soft substrates, Phys. Rev. E 70 (2004) 030601. 

[40] J. Knowles, The finite anti-plane shear field near the tip of a crack for a 
class of incompressible elastic solids. Int. J. Fract. 13 (1977) 611-639. 

[41] O. Kramer, S. Hvidt, and J.D. Ferry, Dynamical mechanical properties 
in: J.M. Mark, B.Erman, and F. Eirich (eds.) Science and Technology 
of Rubber (Academic Press, San Diego, 1994). 

[42] F. John, Plane elastic waves of finite amplitude. Hadamard materials 
and Harmonic materials, Comm. Pure Appl. Math. 19 (1966) 309-341. 

[43] L.D. Landau and E.M. Lifshitz, Theory of Elasticity (Pergamon, New 
York, 1986). 

[44] M. Levinson and I. W. Burgess, A comparison of some simple constitutive 
relations for slightly compressible rubber-like materials. Int. J. Mech. 
Sci. 13 (1971) 563-572. 

[45] G.A. Maugin, Nonlinear Waves in Elastic Crystals (University Press, 
Oxford, 1999). 

[46] F.D. Murnaghan, Finite Deformation of an Elastic Solid (John Wiley: 
New York, 1951). 

[47] A. Norris, Finite amplitude waves in solids. In: M.F. Hamilton and 
D.T. Blackstock (eds.) Nonlinear Acoustics (Academic Press, San Diego, 
1999) pp. 263-277. 

[48] R.W. Ogden, Large deformation isotropic elasticity: on the correlation 
of theory and experimental for compressible rubberlike solids, Proc. R. 
Soc. Lond. A 328 (1972) 567-583. 

31 



[49] R.W. Ogden, G. Saccomandi, and I. Sgura, Fitting hyperelastic models 
to experimental data, Comput. Mech. 34 (2004) 484-502. 

[50] K.J. Parker and R.M. Lerner, Sonoelasticity of organs: shear waves ring 
a bell, J. Ultras. Med. 11 (1992) 387-392. 

[51] E. Pucci and G. Saccomandi, A note on the Gent model for rubber-like 
materials. Rubber Chem. Tech. 75 (2002) 839-851. 

[52] A.V. Porubov, Amplification of Nonlinear Strain Waves in Solids 
(World Scientific, Singapore, 2003). 

[53] P. Rosenau and M.B. Rubin, Motion of a nonlinear string: some exact 
solutions to an old problem, Phys. Rev. A, 31 (1985) 3480-3482. 

[54] P. Rosenau and M.B. Rubin, Some nonlinear three-dimensional motions 
of an elastic string, Physica D 19 (1986) 433-439. 

[55] B.K. Sinha and K.W. Winkler, Formation nonlinear constants from sonic 
measurements at two borehole pressures, Geophys. 64 (1990) 1890-1900. 

[56] J.C. Simo and K.S. Pister, Remarks on rate constitutive equations for 
finite deformations. Computer Meth. Appl. Mech. Eng. 46 (1984) 201- 
215. 

[57] W.H. Steeb and A. Kunick, Lagrange functions of a class of dynamical 
systems with limit-cycles and chaotic behavior, Phys. Rev. A 25 (1982) 
2889-2892. 

[58] C. Truesdell, General and exact theory of waves in finite elastic strain. 
Arch. Rat. Mech. Analysis 8 (1961) 263-296. 

[59] E.A. Zabolotskaya, Sound beams in a nonlinear isotropic solid, Sov. 
Phys. Acoust. 32 (1986) 296-299. 



32 



